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We investigate the distortion of interference fringes in colliding and merging Bose-Einstein Con- 
densates due to inter-atomic interaction effects. We also consider the effect of the distortion on the 
spontaneous formation of vortices. The interference pattern depends crucially on two relevant pa- 
rameters, the relative center-of-mass velocity and peak density of the initial state, which comprises 
a pair of initially well-separated, but coherent, condensate clouds. We identify three qualitatively 
distinct regimes of behavior during the subsequent evolution of the condensates: collision; expan- 
sion; merging. Using a comprehensive set of numerical simulations based on the Gross-Pitaevskii 
equation, we specify the cross-overs between these regimes and identify the system parameters re- 
quired for dynamical instabilities and vortex creation. Interference fringes, which would be planar 
in the absence of interactions, are distorted by the combined effects of interaction and the non- 
uniform density profile of a harmonically trapped condensate. In appropriate parameter regimes, 
the distortion is sufficiently large to produce a net circulation and thereby a dynamical instability 
occurs by which vortices form with their cores on the interference minima. 
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I. INTRODUCTION 

Experiments in which a Bose condensed cloud of ultra- 
cold atoms is divided into two spatially separated com- 
ponents and then allowed to recombine have been widely 
reported [H4lll|. These experiments are of fundamental 
interest, for example, in demonstrating the quantum na- 
ture of the condensate and in investigating decoherence. 
In addition, such processes are central to matter wave 
interferometry using ultra-cold atom condensates, which 
may have many technological applications. In general, 
there are two ways in which initially separated clouds 
can be made to evolve so that they subsequently overlap. 
Firstly, they can be allowed to expand by relaxing the 
confining potential that holds them apart [H, H, 0| • Sec- 
ondly, they can be subjected to external potential gradi- 
ents that cause the clouds to move together and collide 
whilst maintaining their form [Tlj |. I n realistic scenarios, 
a merging process will, to some extent, involve both of 
these mechanisms. The resulting interference pattern will 
be more complex than in either of the idealized cases, par- 
ticularly for systems in which inter-atomic interactions 
are involved. There has been some theoretical work on 
freely expanding condensate clouds in which inter-atomic 
interactions are shown to increase the expansion rate [l2[ 
and create nonuniform interference fringes (l3l . \\a\ . How- 
ever, to date there has been little study of this general 
case or how it compares with either the purely expanding 
or purely colliding scenarios, in particular, the dynamics 
of curved fringes in the merging BECs. 

It is now known that, while the behavior of very 
low density condensates can be well described using el- 
ementary single-particle quantum mechanics, systems in 
which the interactions are stronger show behavior which 
is quantitatively, and frequently qualitatively, different. 



In particular it was shown in Ref.[15| that the interfer- 
ence pattern formed when two countcrpropagating inter- 
acting clouds overlap can give rise to the formation of 
persistent dark solitary waves and, subsequently, the nu- 
cleation of linear arrays of vortex rings. As well as having 
implications for real matter wave interferometers, these 
processes are of intrinsic interest as an experimentally 
controllable route to homogeneous quantum turbulence 
[l6j]. Recently, experiments have observed the sponta- 
neous formation of vortices in a ring trap in which the 
condensate is split into three components and then al- 
lowed to recombine (TtJ- Related theoretical work fo- 
cuses on the role of ramping down the trap potential and 
phase imprinting on the formation of these spontaneous 
vortices [IH, . Although the role of interference in the 
formation of vortices has been pointed out 0, , the 
detailed description of how the interference affects the 
formation of nonlinear excitations has not been shown. 

In this paper, we model experiments in which a Bose- 
condensed cloud is divided adiabatically into two compo- 
nents, for example by using a tailored magnetic poten- 
tial from an atom chip Q , the introduction of an optical 
barrier via a shaped blue-detuned laser Q or by pass- 
ing counterpropagating red detuned laser beams through 
an acousto-optic modulator driven at RF frequencies [3[ . 
We examine how the process of both expansion and col- 
lision affects the interference patterns of two BECs by 
varying the initial peak densities of the clouds and their 
relative center-of-mass (cm.) collision velocities (con- 
trolled by varying the accelerating potential in which 
they move). Our numerical simulations, based on the 
mean field Gross-Pitaevskii equation (GPE), show that 
in regimes where the interactions are significant, interfer- 
ence fringes first develop non-uniform spacings and then 
become curved. In extreme cases, the curvature can be- 
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come sufficient for there to be a net circulation around 
a localized core region leading to the formation of a vor- 
tex ring. Interactions are important at higher densities 
but at lower collision velocities. At higher collisional ve- 
locities the system cannot respond hydrodynamically to 
the high-density fringes before the clouds have passed 
through one another and the fringes have disappeared. 
Through the analysis of a set of simulations and known 
experiments, we argue that a connection between vor- 
tex formation and interference in merging BECs exists. 
Specifically, the interference of the condensates gener- 
ates high-density peaks and extremely low-density val- 
leys. The combined effects of inter-atomic repulsion and 
the non-uniform density profile of the condensates pro- 
duce the additional degrees of freedom in the area of in- 
terference, with the result that the high-density areas 
from the interference peaks encircle the low-density ar- 
eas between them. 

Our general protocol involves preparing each of a pair 
of clouds in the lowest (mean field) energy state of a har- 
monic trap. The prepared clouds are then displaced in 
opposite directions and allowed to evolve when subject to 
a harmonic trap potential, which need not be the same 
as that used to prepare the initial clouds. In general, 
two dynamical processes will occur before the two clouds 
overlap. Firstly the clouds can undergo (ballistic) cm. 
motion in response to the trap potential. Secondly the 
clouds can change shape as they evolve (most notably ex- 
panding or contracting). In general, of course, both pro- 
cesses will occur but we can identify two extreme limits. 
If two copies of a non-interacting Bose gas are prepared 
in the lowest energy state of a given harmonic trap and 
are then placed at symmetrical points away from the cen- 
ter of that trap, the clouds will move ballistically without 
changing volume until they overlap spatially. We refer to 
this as a "colliding process" . Alternatively, two (interact- 
ing) clouds prepared as above can be initially held a fixed 
distance apart and then allowed to evolve after switch- 
ing off the trap potential. In this case, the clouds will 
expand until they overlap. We refer to such a process as 
an "expanding" process. Intermediate cases where nei- 
ther ballistic cm. motion nor expansion dominate will 
be referred to as "merging" processes. 

This paper is organized as follows. Section II formally 
describes the construction of the initial state of two col- 
liding, expanding, and merging condensates, and the nu- 
merical procedure for simulating the dynamics of the two 
condensates. In section III, we derive the interference 
formulae which quantitatively describe the interference 
pattern of two condensates coming together. We also 
identify three distinct dynamical regimes in the cloud 
dynamics. Section IV shows the typical interference pat- 
terns of the three different processes and the formation 
of spontaneous vortices in merging condensates. Also the 
correlation between vortex formation and interference in 
merging process is shown. In section V, we summarize 
general properties in the parameter space. Section VI 
contains our conclusions. 
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FIG. 1: Solid curves: schematic representations of the effec- 
tive potential for trapping two condensates (a) and for the 
merging them (b). Shaded areas in (a), (b) represent the 
atom density profile \&(z, 0, 0)] 2 . 

II. SIMULATIONS 

In the most general case, our simulations have the fol- 
lowing features. We begin with a trap potential of the 
form 

V (r,z) = ^mujl r 2 + ^muj^z 2 , (I) 

where (r, 6, z) are polar coordinates of positions within 
the trap. We choose the form for the function <fio(r, z) 
which minimizes the energy functional 

dZ J drr {^\VM 2 + V (r,z)\ct ) o\ 2 

(2) 

subject to the normalization constraint 

dz / drr \<j) \ 2 rdrdz (3) 

- oo JQ 

where g = 4irH 2 a s /m is the inter-atomic coupling con- 
stant and a s is the s-wave scattering length. We use the 
value a s = 2.9nm and m = 3.82 x 10 _26 kg for a Sodium- 
23 condensate throughout this paper. 

The initial state for the simulations is then obtained 
by setting the order parameter to have the form 

$(r, z, t = 0) = ^ 2(1 ^ Q) {Mr, z-A) + Mr,z + A)) , 

(4) 

where 

/>oo />oo 

Q=—2n dz drrcb (r,z~A)Mr,z + A), (5) 

iV J-oo JO 

representing a pair of clouds displaced in opposite direc- 
tions along the z-axis, as shown in FigfIJa). The param- 
eter Q ensures correct normalization of $ in situations in 
which two atom clouds are not initially fully separated.lt 
is worth pointing out that the representation in terms of 
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a single order parameter, implies that the two clouds 
are fully coherent and it is not the case that there are 
N/2 atoms in each cloud. 

The subsequent evolution of the system is determined 
by numerically solving the Gross-Pitaevskii equation 



ih— = -— v 2 $ 

dt 2m 



V(r,z)<I> + g\<S>\ 2 <I> 



(6) 



assuming that the initial rotational symmetry about the 
z-axis is preserved, using a 2D Crank-Nicholson algo- 
rithm adapted to correctly treat the polar variable r. The 
potential is now 
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V(r,z) 

where in general wn < wmq and ujj_ < lu±.q. 



(7) 



III. GENERAL DESCRIPTION OF TWO 
MERGING CONDENSATE CLOUDS 

In this section, we analyze the interference pattern gen- 
erated by the merger of two condensates and identify 
three types of behavior and their corresponding parame- 
ter regimes. We refer to the two extreme cases as colli- 
sion and expansion, and to the generalized intermediate 
case as merging. While the interference fringes described 
by the (linear) Schrodingcr equation arc planar, the in- 
teractions included in the GPE via the non-linear term, 
o|(j)| 2 (j) ; give rise to phenomena such as curved fringes 
The repulsive interaction depends on the density 
distribution in the interference region of the two conden- 
sate clouds. Consequently, there is a position-dependent 
rate of expansion, resulting in different fringe spacings at 
different positions. 

Our purpose is to explore the relationship be- 
tween nonlinear interactions and the expansion rate 
of condensate peaks, which we quantify using rj = 
r o °° z\<f>(z,t)\ 2 dz/ f °° z\<S>(z,0)\ 2 dz. The system of inter- 
est is a series of condensate fringes with similar shapes 
but different densities. The presence of interactions 
should cause the fringes to expand at different rates, and 
hence their spacings to change. In order to approximately 
simulate the expansion of individual fringes, we consider 
the ID model of a freely expanding single condensate 
initially held in a harmonic trap. To consider different 
fringe densities, the numbers of atoms in the model are 
changed, and then the trap frequency adjusted to create 
a similar density profile. Inspired by a 3D model of a clas- 
sical gas (20[, we introduce a scaling factor r/(t) for our ID 
case, where the condensate width R(t) = r](t)R(0). The 
analytic solution of the expansion rate can be derived by 
adapting the method used in [2(| to the ID case. In the 
Thomas-Fermi approximation, n(t) satisfies the dynamic 
equation [2(| 
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where u> is the trap frequency. By integrating Eq. 
find 



2 (Vv(v - 1) + log 1 07 + Vv^ T|) = 



For small t, the solution for Eq.© is 

If t is large, rj » 1, then we obtain 
rj = \ 2ojt. 



(9) 
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Eqs. fTUl) and (TTTj) indicate that a higher peak density, 
corresponding to larger w, causes a higher rate of expan- 
sion in ID freely expanding clouds. Fig 12] shows that our 
simulations are in basic agreement with analytic predic- 
tion that a higher density leads to a larger expansion 
rate. Although we fit the numerical data well at small 
and large t respectively using generic parabola and lin- 
ear formulae, the fitting formulae are not the same to 
Eq. ([T0|) and Eq . (fTT]) . This originates that at short time, 
the expansion of the cloud is speeded up due to the inter- 
atomic interaction, and then the Thomas-Fermi shape of 
the density distribution is broken down. 
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FIG. 2: The time-dependent expansion rate of the ID freely 
expanding condensate, with gm = 1.3a s , from our simulation 
(curves from top to bottom have N = 100, 400, 600, 800); 
Insert: short-time behaviors. 

Furthermore, we take into account the effect of the 
repulsive interaction on interference fringes by introduc- 
ing a variable a, which is a modified coefficient depend- 
ing monotonically on the expansion rate of interference 
fringes. Qualitative insight into the properties of the in- 
terference of two condensate waves can be gained by us- 
ing a Gaussian ansatz, 



,t) = C 




-(l-i X (t))z 2 /2^ 2 (t) 



(12) 



for the wave function, where C = (ttIq 



,-V4 



x(t) 



ht/ml 2 , + at I ml and a 2 (t) = Z§ (l + X 2 (t)). 

Now, we arrange a pair of wavepackets initially at po- 



sitions ±L/2, $(z,0) 
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(0(z-f,O) + 0( Z + §,())). 
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Then they move toward each other with the same speed 
v = hk /m, described by 



$(z,t) = 



1 . , L 
-i=(<Kz--,t)e 



ik ' z + <t>(z+-,t)e ik '*), (13) 



where the wavevector k is related to the cm. speed of 
atom cloud, depending on the variation of trap potential. 
Inserting Eq. ([T2")) into Eq. (fT5|) . we deduce the expression 
of |<I>(z,i)| 2 , and the cosine term in |$(z,i)| 2 shows the 
effective wavelength for the interference fringes, 
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Assuming that <7i£>|$| 2 is zero (i.e., single-particle pic- 
ture) and the shape of the clouds is frozen, the time- 
dependent terms in Eq. fH)) vanish and thus A = p- = 
— , analogous to classical interference. 

Setting, k! = so that the centers of the initial states 
don't move and (?id|$| 2 = 0, means that the two wave 
packets expand freely and the interference fringe spacing 
then is 



h 2 t 2 

A = 2n(l^m-\ )/HtL 



(15) 



We refer to this regime as expanding. In general, the first 
term in the bracket in Eq. (|15|) can be neglected with re- 
spect to the second term at large t, resulting in the fringe 
spacing A = ht/mL, which has been used in experiment 
|l| to explain the interference of two expanding conden- 
sates. 

Except for these limiting cases, all other processes 
of two condensates coming together are in the merging 
regime. Obviously the inter-atomic interaction affects 
the interference pattern in merging process. For exam- 
ple, if the role of the nonlinear interaction is involved in 
the collision, Eq.fTTJJ becomes 
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h 2 



(16) 



where the larger inter-atomic interaction, the larger a, 
leading to the larger A. 



IV. DISTORTION OF FRINGES IN TWO 
INTERFERING CONDENSATES 

In this section, we illustrate numerically "collision", 
"merging" and "expanding" behaviors of two BECs. In 
particular, we focus on the role of the repulsive inter- 
atomic interactions in the distortion of fringes and the 
resulting spontaneous vortex formation. 



A. System I: Free Expansion of Well Separated 
Clouds 

We set the initial separation of the two condensates, 
A = 4.8/11 with Z|| = yJh/muj\\_Q, where each condensate 



wavefunction </> is the groundstate solution of Eq.© with 
trap frequencies cjiiq = 27r x 180Hz and Wj_,o = 27r x 
120Hz. Then, the system is allowed to evolve freely by 
setting u)\\ = lo±_ = 0. In the absence of interactions, the 
order parameter has the form $(r, z — A, i) + <J>(r, z + A, t) 
where 21 
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ttHJ ^ W V ' (17) 

X e -("||(*)-^||(*))z 2 /2Zy e -(ai-*/3i)r 2 /2il 



with a>j = 1/(1 +w 2 t 2 ) and /3j = ujjt/ (l + w|i 2 ) . 
Hence at arbitrary time t we have 

n(r,z,t) = Mr,z,t)\ 2 

- 2a J- (*)"!! ( f ) p -a I (t)r 2 /l 2 , -a„( Z 2 +A 2 )/i 2 



?2l^i + c„s « (18) 



x cosh 



From the cosine term in Eq . p8|) . we deduce that the 
length scale characterizing the interference fringes is 



A(t) 



nlf { 1+ujft 2 
A wtit 



(19) 



which, for t > 2n/(J\\, behaves as A(t) ~ (ttI 2 / A)ui\\t 
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FIG. 3: Gray-scale plots of atom density (black=high) in the 
z-r plane (axes inset) for double condensates, evolving with 
the process of expanding at t = 5.81ms(a) and 10.21ms(b). 
and 1 label respectively zeroth- and first-order fringes. 

In FigUJa), our simulations show that, at very low 
densities (N = 3 x 10 4 ), the growth of the fringe spac- 
ing is the same for the zeroth- and first-order fringes 
(depicted in FigJ^b)) and linear in t, which matches 
the non- interacting expression well 0, [ljj . At a slightly 
higher density (N — 4 x 10 5 ), the interference dynam- 
ics exhibits three regimes, as indicated in Figf4] (b). In 
regime I (t < 2.25ms), the growth of the fringe spacings 
is identical for the zeroth- and first-order fringes due to 
the low densities in the initial merging region. Note that 
this is distinct from the non-interacting behavior because 
the nonlinear interaction increases the expansion rate of 
the clouds [l2|. In regime II (2.25ms< t < 9.75ms), the 
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spacing of the zeroth-order fringe grows faster than that 
of the first-order fringe, indicating that the ideal fringe 
pattern is distorted in the z-direction. Specifically, ow- 
ing to its higher density, the central fringe is thickened 
relative to its neighbors. In regime III (t > 9.75ms), 
the growth rate of the zeroth- and first-order fringe spac- 
ings reduces gradually, approaching the non-interacting 
limit. For enough long times, the difference in the zeroth- 
and first-order fringe spacings becomes negligible and the 
fringes are uniform. 
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FIG. 4: The zeroth-order (solid line) and the first-order 
(dashed line) fringe spacings of two freely expanding conden- 
sates versus time t with N = 3 x 10 4 (a) and N = 4 x 10 5 (b). 
The time evolution of the peak densities of the zeroth- (dot- 
dashed curve with circles) and first-order (dot-dashed curve 
with crosses) fringes (c). Dotted line in (a) labels the time, 
t = 3.3ms, at which the zeroth- and first-order fringes appear 
and dotted lines in (b) label the times that separate the three 
interference regimes. 



B. System II: Collision of Well Separated Clouds 

In this set of simulations, the initial state is prepared 
as before for a range of N. However, the trap remains 
present during the dynamical process, as depicted in 
FigOJb), so that the two components are accelerated 
towards one another. For low-density condensates, the 
system behaves, at least over one period of the longi- 
tudinal trap, like an ideal Bose gas. Since wj_ = wj_.o 
and W|| = wy 0) the initial state is a superposition of two 



coherent states of the oscillator potential, so that 
$( r , Zjt ) = _L Ce-^' 2l ^(e- iK ^e-^- z ^' 2l l 

+ e iK(t)z e -(z+Z(t)) 2 /2l\^ 



(20) 



where Z(t) = A cos (^J\\t) and K(t) = — sin (wyi). When 

the clouds are maximally overlapping (t = 7r/2a»y), the 
fringe spacing is given by, 



2tt 27r/ff 



Kit) 



(21) 



(a) 




(b) 


• 


• 


III 



(e) 



2 - 





(g) 




5pm z 



(d) 




m 





FIG. 5: Gray-scale plots of atom density (black=high) in the 
z-r plane (axes inset) for double condensates, evolving with 
the process of colliding at t — 0ms (a), 2.01ms (b), and merg- 
ing at t = 1.21ms (c), 2.01ms (d), 2.41ms (e), 3.61ms (f), 
4.11ms (g), 4.41ms (h); The phase of the BEC wave function 
within the dashed box in (h) is shown in the gray-scale plot 
(black=27r, white=0) to the bottom of the density profile. 
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In our simulations, the initial density distribution is 
shown in FigJSJa). Subsequently the two clouds approach 
and collide at speed v given approximately by 

\riv 2 = A 2 , (22) 

which gives V = wyA and A = h/mv = 27r/ 2 /A. We 
find that the collision velocity is V — 6.08/||W|| = 1.07 x 
lO^ms" 1 and A = 1.036Z|| = 1.62 x l(T 6 m. 
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FIG. 6: The peak density of the zeroth- (dashed curve) and 
first-order fringes (solid curve), and the zeroth- (solid curve 
with crosses) and first-order fringe width (dashed curve with 
circles) versus the initial peak density of the condensate, no, 
in "colliding" process with fixed g. 

Fig|6] shows the effect of the nonlinear term on the 
fringe pattern of the colliding condensates with differ- 
ent atom numbers. For an initial peak density uq < 
1.52 x 10 14 cm -3 , the peak density of the zeroth- and 
first-order fringes increases linearly with atom number, 
as expected for interference in a system governed by a 
linear wave equation. In this regime, for given initial trap 
frequencies and separation of the clouds, the fringe spac- 
ings are independent of the initial density distribution. 
However, when no > 1.52 x 10 14 cm -3 , the peak densities 
of the two fringes depend sublincarly on no. The gap be- 
tween the two fringes increases with no, with the result 
that the differences between fringe spacings increase as 
well. When this difference in the high-density region be- 
comes sufficiently large, the fringes become significantly 
curved leading, ultimately to vortex formation. 

Furthermore, our theoretical formula Eq.lfTHl) shows 
a proper explanation for the nonuniform fringes. The 
higher density in the zeroth-ordcr fringe, characterized 
by larger a, leads to a larger fringe spacing than that of 
the first-order fringe. 

C. System III: Merging of Well Separated Clouds 

Here, we choose the initial state to have no = 4.6 x 
10 21 m -3 and A = 4.33 x 10 _6 m. The initial wave func- 
tion 4>o is obtained by numerically minimizing the mean 



field energy Eq.Q with the trap frequencies, wm o = 
2tt x 800Hz and uj ±fi = 2tt x 533Hz. Then at time t'= 0, 
we change the trap frequencies wm = 2ir x 180Hz and 
uj± = 2ir x 120Hz. Since the trap used for the evolution 
of the two clouds is much less confining than the trap 
used in preparing them, each cloud expands due to inter- 
nal pressure as well as undergoing bulk motion. A here is 
smaller than that in system II and hence the relative ve- 
locity of the two clouds at maximal overlap, v = 2.8Ziia;i|, 
is slightly smaller than in system II. 

FigEJc)-(h) show the time evolution of the density pro- 
file. The process of both expansion and collision leads 
to a distorted fringe pattern with the characteristic con- 
figuration that there are "thicker" fringes in the center 
of the cloud and a larger fringe spacing, towards the 
edge of the cloud where interaction effects are negligi- 
ble and the fringe spacing is closer to its non-interacting 
value (FigJSfc)). This curved fringe pattern is similar to 
that observed in experiment [l[ where condensates with 
short separation tend to form large variations in den- 
sity in the merging region. After the clouds have com- 
pletely merged, the fringe spacings become smaller, and 
the distinction between the center and edge part as well 
as the zeroth-order and first-order peak, reduces corre- 
spondingly (FigJSfd)). 

The complex interference pattern can also be inter- 
preted by Eq. (fT4"|) . At the beginning of the merging 
(t < 2.1ms), the expanding speed of an atom cloud v e 
is larger than its cm. speed, v m , due to very large no, 
and thus the fringe spacings arc dominated by the term 
at/h in Eq. (fl~4"| . It implies that the larger peak density of 
the fringes, the larger a, resulting in larger A. Therefore, 
the fringe spacings approaching z = 0, are slightly larger 
than those at the edge, originated from the peak density 
of the zeroth-order fringe larger than those of high-order 
fringes. With time evolution in the merging process, the 
increasing k and ht/2m in Eq. (|14[) gradually dominate, 
resulting in smaller fringe spacings and less variation of 
the fringe spacings between the middle and edge of the 
clouds. 

As shown in FigJSJd) and (e), the strong interactions in 
the center of the clouds leads to a net radial flow of atoms 
in the high density fringes, triggering a sound wave in the 
central peak to propagate radially from position 1 to 2. 
Our calculation shows that v r w £1.2/^1,2 = ll.lmms -1 , 
larger than 2v m (« 9.8mms _1 ). This rapid outflow leads 
to complete depiction of the center of trap after the two 
clouds have passed through one another, resulting in the 
formation of a dark soliton (e). When the two clouds 
recollide at time 2t, both are recollapsing radially lead- 
ing to an enhanced density in the bright fringes causing 
an even greater longitudinal expansion and greater den- 
sity gradients in the radial direction (f). The curvature 
of the fringes is thus even greater than in the first col- 
lision, (g), and is sufficient to generate a net circulation 
around localized regions that are fully depleted: vortices 
are generated, as seen in (h). 

In general, due to strong inter-atomic interaction and 
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the competition between v m and v r , additional degrees 
of freedom are generated, with the result that the high- 
density areas from the interference peaks encircle the low- 
density areas between them (FigJSJg)). The positions 
of the vortex rings correlate naturally with the previous 
positions of interference valleys. Naturally, the formation 
of vortices from quantum reflection of high-density and 
low-velocity BEC [lj| is analogous with our results. 



V. GENERAL FEATURES IN THE 
PARAMETER SPACE 

To explore the interference properties and characterize 
the spontaneous formation of vortices over the param- 
eter space, we have performed a comprehensive set of 
numerical simulations to obtain FigJT] Each data point 
in FigJT] (a) and (b) is derived from a set of at least 12 dy- 
namic calculations. We identify the three types of process 
through their interference characteristics. In an expan- 
sion process, the interference fringes are time-dependent 
and uniform. In the colliding process, the fringes are 
static and uniform while the clouds are fully overlap- 
ping. The interference pattern in a merging process is 
time-dependent and nonuniform. 

FigUJ (a) demonstrates that starting in the freely ex- 
panding regime, moving to higher initial densities or 
smaller initial separations will produce the more gen- 
eral merging behavior. This is exhibited by a loss of 
homogeneity in the fringe spacing and, in turn, fringe 
curvature. The upwards curve of the boundary between 
expanding and merging behavior is similar to that for 
the boundary between colliding and merging behavior 
shown in (FigUJ (b)). The rigorous explanation of this 
trend might be derived from the quantitative analysis of 
the effect of nonlinear interaction on the fringe spacing, 
which is beyond the range of this paper. 

Based on the analysis of the dynamics over a large 
portion of the parameter space, FigJT] (c) summarizes the 
formation of vortices in merging BECs. Above a critical 
initial density, n c , vortices will form during the return 
overlap of the two clouds within a definite range of initial 
cm. velocities. If the velocity is too low then damping 
of the motion during the first overlap effectively arrests 
the dynamics. In addition, the low velocity leads to a 
large fringe spacing so that, even towards the edges of 
the cloud, the curvature is never sufficient to cause a net 
circulation. If the collision velocity is too high, the hydro- 
dynamic response of the cloud is too slow and the clouds 
have separated before any significant distortion occurs 
[22^ (inserted plot in FigUJ ( c ))- Our simulations show 
that the curvature of the interference fringes is of impor- 
tance in the formation of vortices and that the number 
of fringes within the cloud determine the number of vor- 
tices. This is shown in the inserted plots of Fig.(|7j) (c) 
where the five fringes from the short initial separation 
result in the formation of two vortex rings while the nine 
fringes from the large initial separation cause six vortex 
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FIG. 7: In (a) and (b) the cross-overs between different qual- 
itative behaviors are illustrated. The generic "merging" pro- 
cess is identified by a fringe non-uniformity of at least 0.4/y 
between zeroth- and first-order fringes in case (a) and at least 
10% in case (b). In (a) a pair of identical clouds are placed 
a distance A apart and allowed to expand in the absence of 
any trap potential. In (b) the clouds evolve in a single trap 
potential. In (c) the locus of points in the parameter space 
above which vortices are produced is shown schematically for 
a generic process. Insert: three density configurations for the 
same clouds with different separations in the merging and 
overlapping process with their parameter points arrowed. 



rings. The connection between the distortion of fringes 
and the formation of vortices can explain the experiments 
in which more vortices are generated in the faster merg- 
ing, interfering region |(|. The faster merging produces 
more interference fringes, creating the possibility for the 
formation of more vortices in the region. The diagram 
in Fig.© is based on the specific material parameters 
described in the previous sections. Simulations using dif- 
ferent parameters give results that are similar but with 
the cross-over lines shifted. In general, for the smaller 
values of the coupling constant, g, and with the larger 
numbers of atoms, it is easier to observe the crossover 
from colliding process to merging process. 

Finally, we address the effect of the halo of scat- 
tered atoms, produced by the counter-propagating con- 
densates, on the interference pattern. The phenomena 
that two BECs collide at a sufficiently high velocity, a 
halo of elastically scattered atoms is produced, have been 
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demonstrated by both experiments [23l-l25j and theory 
prj [27j . In the expanding and merging processes of our 
system, the scattered atoms might not have significant 
effect on the interference patterns because the rapid de- 
crease of the densities corresponding to the fast expansion 
of the clouds would reduce the rate of atom-atom colli- 
sion in merging areas. For the colliding process, the vis- 
ibility of the interference can be reduced and the fringes 
can be distorted due to the scattered atoms. However, 
these effects can be suppressed or avoided if the density 
and velocity of clouds are low, or the collidingprocess is 
controlled by properly adjusting the barriers [a. ITU . 

VI. SUMMARY 

In summary, we have investigated vortex formation 
and interference instability of two merging condensates 



and identified three distinct regimes in situations rele- 
vant to atom interferometry. Our simulations show an 
explanation about recent experimental work Q, where 
the faster merging three BECs creates more vortices. The 
regularity of vortex formation and interference might al- 
low for the design of experiments to study vortex creation 
and the dynamics of regular, linear, vortex arrays. Our 
latest calculations show clearly that one could control the 
density of vortices generated by tuning the fringe spacing. 
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